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Abstract 

We derive basic equations of electromagnetic fields in fractal media 
which are specified by three indepedent fractal dimensions on in the re- 
spective directions Xi (i = 1, 2, 3) of the Cartesian space in which the frac- 
tal is embedded. To grasp the generally anisotropic structure of a fractal, 
we employ the product measure, so that the global forms of governing 
equations may be cast in forms involving conventional (integer-order) in- 
tegrals, while the local forms are expressed through partial differential 
equations with derivatives of integer order but containing coefficients in- 
volving the cti's. First, a formulation based on product measures is shown 
to satisfy the four basic identities of vector calculus. This allows a gen- 
eralization of the Green-Gauss and Stokes theorems as well as the charge 
conservation equation on anisotropic fractals. Then, pursuing the concep- 
tual approach, we derive the Faraday and Ampere laws for such fractal 
media, which, along with two auxiliary null-divergence conditions, effec- 
tively give the modified Maxwell equations. Proceeding on a separate 
track, we employ a variational principle for electromagnetic fields, appro- 
priately adapted to fractal media, to independently derive the same forms 
of these two laws. It is next found that the parabolic (for a conducting 
medium) and the hyperbolic (for a dielectric medium) equations involve 
modified gradient operators, while the Poynting vector has the same form 
as in the non-fractal case. Finally, Maxwell's electromagnetic stress tensor 
is reformulated for fractal systems. In all the cases, the derived equations 
for fractal media depend explicitly on fractal dimensions and reduce to 
conventional forms for continuous media with Euclidean geometries upon 
setting the dimensions to integers. 



1 Introduction 

In a recent article (Li & Ostoja-Starzewski, 2009) we formulated a continuum 
mechanics of anisotropic fractal elastic media using two independent approaches: 
a mechanical (Newtonian) one and a variational (Lagrangian-Hamiltonian) one. 
Proceeding separately on both paths, we derived the identical equations of 
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wave motion in one-, two- and three-dimensional media, and this provided a 
verification of the consistency of the product measure capable of grasping the 
anisotropic fractal media as well as the verification of the final equations. The 
question arises whether one can use the same strategy to derive the field equa- 
tions of electromagnetism in fractal media, be they isotropic or anisotropic. 
Various studies generalizing the classical, quantum or relativistic field theories 
to fractal spaces (Stillinger, 1977; Svozil, 1987; Palmer & Stavrinou, 2004; Not- 
tale, 2010) as well as the electromagnetism to fractal media appeared over the 
past two decades (Tarasov, 2006, 2010). They have been formulated either 
through a conceptual path — in the case of electromagnetism, by developing 
the Gauss, Faraday and Ampere laws, so as to obtain the Maxwell equations — 
or through variational principles. However, a derivation obtained on one path 
(say, conceptual) has never been verified through another, separate derivation 
(respectively, variational) . 

In the following we proceed under the premise that the conceptual derivation 
of Maxwell's equations for anisotropic fractal media should yield the same results 
as the derivation based on a variational principle, both derivations being based 
on the same dimensional regularization. To grasp the generally anisotropic 
structure of a fractal, we employ the recently introduced product measure, so 
that the global forms of governing equations may be cast in forms involving con- 
ventional (integer-order) integrals, while the local forms are expressed through 
partial differential equations with derivatives of integer order but containing 
coefficients involving fractal dimensions on in three respective directions Xi, 
i = 1,2,3. 

We first show that the formulation based on product measures satisfies the 
four basic identities of vector calculus. This readily leads to a generalization 
of the Green-Gauss and Stokes theorems involving the fractal gradient and curl 
operators, and hence to the charge conservation equation on anisotropic fractals. 
Then, pursuing the conceptual approach, we derive the Faraday and Ampere 
laws for such fractal media, which, along with two auxiliary null-divergence 
conditions, effectively give the modified Maxwell equations. 

Proceeding on a separate track, we employ a variational principle for electro- 
magnetic fields, appropriately adapted to fractal media, to independently derive 
the same forms of these two laws. We next examine whether the parabolic (for 
a conducting medium) and the hyperbolic (for a dielectric medium) equations 
involve fractal gradient operators, and also whether the Poynting vector has the 
same form as in the non-fractal case. Finally, Maxwell's electromagnetic stress 
tensor is reformulated for fractal systems. In all the cases, the derived equations 
for fractal media depend explicitly on fractal dimensions and reduce to conven- 
tional forms for continuous media with Euclidean geometries upon setting the 
dimensions to integers. In order to make the presentation clear, in some places 
we give the tensorial relations in the symbolic as well as the index notations. 
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2 Background 



2.1 Anisotropic Fractal Distributions 

The distribution of charges in an electrically and magnetically conducting isotropic 
fractal medium embedded in the Euclidea space E 3 obeys a power law 

Q(R) ~ R D , D<3, (1) 

where D is the fractal dimension and R is the spatial resolution. However, in 
an anisotropic fractal medium (generally denoted by W) we write 

Q(x u x 2l x 3 )~x 1 a ix 2 a2 x 3 a3 , (2) 

where each a k plays the role of a fractal dimension along each axis Xk {i — 1,2,3) 
of the Cartesian space in which the fractal is embedded. In other directions, 
the fractal dimension is not necessarily the sum of projected fractal dimensions, 
however, as noted by Falconer (2003), "Many fractals encountered in practice 
are not actually products, but are product-like. ." Hence, we expect the equality 
between the fractal dimension D of the total charge and ct\ + a 2 + 0:3 to hold 
for fractals encountered in practice. When D — > 3, with each a.i — > 1, the 
conventional concept of charge distribution is recovered. 

By analogy to the formulation of continuum mechanics on anisotropic frac- 
tals (Li & Ostoja-Starzewski, 2009), we now express (2.1) through a product 
measure 

Q(W)= p(x 1 ,x 2 ,x 3 )dl ai (x 1 )dl a2 (x 2 )dl a3 (x 3 ), (3) 
Jw 

where p is the density of distribution. In (2.3), the length measure in each 
coordinate is given in terms of the transformation coefficients c[ k ^ by 

dl ak (x k ) = c[ k \a k ,Xk)dx k , A; = 1,2, 3 (no sum), (4) 

where l k is the total length (integral interval) along x k and l k o is the character- 
istic length in the given direction, like the mean pore size. The relation (2.4) 
implies that the infinitesimal fractal volume element, oIVd, is related to the cubic 
volume element, dV 3 — dx\dx 2 dx 3 by 



dVo = dl ai (xi)dl a2 (x 2 )dl a3 (x 3 ) — c^ 1 c^ 1 dx\dx 2 dx 3 — c 3 dV 3 , 
with c 3 = c^c^c^. 



(5) 



Similarly, the infinitesimal fractal surface element, dSd, is related to the planar 

(k) 

surface element, dS d — dxidxj with the normal vector along x k , according to 

dSf ] =c c dS 2 , 

.,, (fe) (») (j) 1 (fe) . , . • • / , w 

With C 2 ' =c\'c\ J > =C 3 / C \>, l^J, l,J j£k. 
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The sum = cei+aj, i 7^ j, i, j 7^ k, is the fractal dimension of the surface 
along the diagonals | x 2 1 — I X j I in S d ' . Again, this equality is not necessarily true 
elsewhere, but is expected to hold for fractals encountered in practice (Falconer, 
2003). 

The transformation coefficients showing up in (2.4) can be represented 
in terms of the modified Riemann-Liouville fractional integral of Jumarie (2005, 
2009), thus 



(fe) 



Ik %k 

ho 



ctk — 1 



k = 1,2, 3, (no sum), 



(7) 



although further developments will not depend on this explicit form. 

2.2 Vector calculus on anisotropic fractals 

Following (Li & Ostoja-Starzewski, 2009), we will extensively employ the fractal 
derivative (fractal gradient) operator (grad) 



V D cj> = e fe V£V or V£\ 



»(*) dx k 



(no sum on k). 



(8) 



where are base vectors. Hence, the fractal divergence of a vector field 

*f-v».f or v?/ s = -Lg. 



(9) 



Note that this leads to a fractal curl operator of a vector field 

S k ) dx k ' 



curlf = V D x f or x /, = e jkl 



D 



(10) 



We observe that the four fundamental identities of the conventional vector 
calculus carry over to these new operators: 

(i) The divergence of the curl of a vector field f : 



div.curl / = Vfxf= ' 



S k ) dx k 



(ii) The curl of the gradient of a scalar field 0: 



curl x (grad 4>) = Vf x (Vf x<p) = e ijk 



_1 d_ 



1 1 dfj 
£jki c U) C W dXj dx k ~ U - 

(11) 
1 1 df t 



c[ k) dx k 



Gjki 



(12) 

In both cases above we can pull 1/c^ in front of the gradient because the 
coefficient is independent of Xj. 



0. 
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(iii) The divergence of the gradient of a scalar held 
the fractal gradient as 



<f> is written in terms of 



?D V7-D, 



1 d 


1 d(j) 


1 


d<p,j 






^dx~_ 


c i 
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(13) 



which gives an explicit form of the fractal Laplacian. 
(iv) The curl of the curl operating on a vector held f : 



curl x (curlf) = V^ 5 x (V, D x f ) = e f 



1 d_ 

S r ) dx r 



( : jki 



^ dxk 



_1 d_ 



J_djj_ 

L c[ h) 9x k 



(SkpSir - 5kr 



1 d_ 

» dx r 



J_df^ 

S k ) dx k 



1 d 


' 1 df r ' 


1 d 


' 1 df p 


c ( r ) dx r 


_ c [p) dx p _ 


C W dx r 


c[ r) dx r 



v? (vf • f) - V?V? f 

(14) 

As further background, we now extend two integral theorems to generally 
anisotropic fractals. 

2.3 Stokes and Green-Gauss theorems for anisotropic frac- 
tals 

Stokes theorem. We begin with J A (V D x f) • ndSd and proceed in the index 
notation: 

/ n k Vf f t dS d = ( n k e kjl -7-^fi,j dS d = [ n k e kji -7-r fi,j c 2 k ' ' dS 2 

JdW JdW c\' JdW Ci 

= / n k e k jifi,j c^dS 2 = / n k e k ji (fiC^) ,j dS 2 

JdW JdW V 7 

(15) 

In the second equality above we employed the dimensional regularization, fol- 
lowed by taking note of (2.6), including the fact that cf 1 is independent of xj. 
Next, by the Stokes theorem in the Euclidean space E 3 , we have 

/ n k e kji (fiC^) ,j dS 2 = \ f l c i ^dl i = \ f l dl lD (16) 
JdW v ' JdW J A 

where, again, we used (2.4). Thus, we arrive at the Stokes theorem for fractals 

n- curlf dS d = [ f-dl ai (17) 
Ji 



Formally, the above agrees with Tarasov's (2006) result in the case of isotropy, 
although our curl operator also grasps anisotropy. 
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Green-Gauss theorem. We begin with J A f <g> ndS d and proceed in the 
index notation: 



/ f t n k dS d = / fin k c ( 2 k) dS 2 = / (fiC 2 k) ), k dV 2 
Jaw Jaw Jw v 7 

/ f/»4 fe) ) 'fc c 3 = / fi,kC 2 k) Cs l dV D = ( fi,k-^xdV D , 
Jw v 7 iw iw c ; 



whereby in the second equality we employed the dimensional regularization, 
followed by taking note of (2.5), then followed by employing the Green-Gauss 
theorem in E , and in turn followed by noting that c 2 i s independent of x k . 
Thus, we arrive at the Green-Gauss theorem for fractals as 

/ t-ndS d = [ V D t dV D (19) 

JdW Jw 

The same comment as that following (2.7) applies here. 

2.4 Charge conservation on anisotropic fractals 

Let J be the current density, n be the direction normal to the surface, and rj be 
the charge density on the fractal W. Then we have 

/ J • n dS d = - [ T] dV D (20) 

JdW Jw 

On account of (2.18), we obtain 

global form f V D -3dV D = -^-[ rj dV D 



d 

local form V D ■ J = — — n 

dt 1 



(21) 



The fact that S7 D automatically appears in the above will be exploited below. 

Focusing on a linear electromagnetic response, with er being the electric 
conductivity tensor, Ohm's law for an anisotropic medium reads 

J = er • E or J % = atjEj (22) 

the isotropic form following for cr^ = S^a. The fact that the constitutive 
form above carries over directly from the non-fractal case is motivated by the 
analogy to elastic media, where Hookc's law is unchanged when going from non- 
fractal to fractal media (Li & Ostoja-Starzewski, 2009); that result ensured the 
consistency of the Newtonian and Hamiltonian approaches to the derivation of 
governing equations - recall the beginning of the Introduction. 
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3 Formulation of Maxwell's equations via Stokes 
and Green- Gauss theorems 

3.1 Faraday's law 

First, note that the Faraday law holds for a fractal W embedded in E 3 



d_ 
~dt 



J^B- ndA d = - jf E ■ dl D (1) 

On account of the Stokes theorem for fractals (2.17), we obtain 

d (' d f d f 1 

di J A BkUkdAd = dt J A ek ^f E * nkdAd = ~~fo J e k 3 t-(jj E il3 n k dA d (2) 

which, by localization leads, to 



d d 1 



or 



0= 4-B + V D x E 

dt 



(3) 



3.2 Ampere's law 

First, note that the Ampere law holds for a fractal W embedded in E 3 

/ C • ndA d = [n-dl D = In- ( V D x H)dA d (4) 

J A Jl Jl 

Noting 



we obtain 



C = J + f (5, 



•I • ^ V" II . ((i) 



which, again by focusing on the integrand, yields 

J + —- V"° x H = (7) 

dt 

On account of the speed of light being related to the dielectric (eo) and magnetic 
(^ ) constants by the well known relation c = * ne above is equivalent 

to 

1 j + <*L _ c 2 V ^ x B = (8) 

e at 
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The equations (3.7) and (3.8) are subject to the constraints 



Gauss law for magnetism V D • B = 
Gauss law V D • E = 



D v n (9) 



where the presence of the fractal (rather than the classical) divergence operator 
\7 D has clearly been suggested by (2.21). 

4 Derivation from variational principle 

The variational principle for the Maxwell equations in conventional (non-fractal) 
setting E 3 



V-E = 0, VB = 0, 

$B <9E 1 

— +VxE = 0, — -c 2 V x B + -J = 0, 

at at e 



(1) 



is (Seliger & Whitham, 1968) 



CdV 3 dt = 0, 

with ( 2 ) 



eo 



11/ dE- 1 
-c 2 B t Bi - -EiEi + Ai \ -q^ - c 2 SijkBk,j + — Ji ) + X E %, 



Here the vector potential A and the scalar potential \ are the Lagrange mul- 
tipliers, while the integrand of (4.2) is slightly extended to include the current 
density J. By varying (4.2)i with respect to E, B, A, and x, gives, respectively, 

8E : E = -dA/dt - Vx 

SB : B = V x A 

1 (3) 

5A : dE/dt — c V x B H J = 

eo 

5x : V • E = 

It is clear that (4.3)1,2 satisfy the first pair of Maxwell's equations (4.1) identi- 
cally. Henceforth, our goal is to obtain the Maxwell equations (4.1)i-(4.1) 4 for 
an anisotropic fractal, by using the same type of a variational approach. 

For a fractal W embedded in E 3 , the variational principle for electromagnetic 
fields is written as 

r t 2 r 

£dV D dt = (4) 



'77 

Jti Jw 
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with C given by (4.2)2- 

Varying (4.4) with respect to E, we obtain 



E 



= -dA/dt- V D ( XCl ) or Ei = -dAi/dt— (xc[ l) ) » = Vf (a^) 

c i 

(5) 

while varying it with respect to B, we obtain the Gauss law for magnetism on 
fractals 

B = V D x(A Cl ) or B k = e iJk -^(A i <% ) ), j = V?(Ai($ i) ) (6) 



In both cases, the fractal gradient operator, V 15 , shows up due to the presence 
of dV D in (4.4) i; recall (2.5) i. 

The next step is to verify whether substituting (4.5) into (3.3) would give a 
zero vector. Thus, 



9 1 1 

= -Qj_B k + ehji-^Ei,j = dB k /dt + e kjl — 



.(<) 



|K»)-(x4-) 



dB k /dt + e k]l — 



d_ 

dt 



-(Ac (i) ) -(yc {i) ) 



d_ 

dt 



Bu-Vf x(^c«)]-Vf x (*<£>)„ 



(7) 

The first term in the square brackets in the last line of the above agrees with 
(4.6), while the last term vanishes, providing the fractal is isotropic (c^ = 
i =/= j), in which case (4.7) is satisfied identically. Thus, for an anisotropic 
fractal, there appears a source/disturbance — an observation consistent with 
Tarasov's formulation, though his explicit form differs from ours. 
Also, substituting (4.6) into (3.9) i yields 



° = e ^~W) 4) ( Aic([j) ) = v ^ v ? x ( Atc({ 



0) 



c i c l 



(8) 



which vanishes providing the fractal is isotropic (c^ = i ^ j). 
Varying (4.4) with respect to A, we obtain 

d~E/dt — c 2 V D x B + — J = or dEjdt-c^f (B k c[ j) ) + — J t = (9) 

which agrees with (3.8) perfectly. 

Varying (4.4) with respect to x, we obtain the Gauss law 



V D -E = or V?E i = -^E i , i = 



(10) 
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which perfectly agrees with (3.9)2- Thus, we have derived by an entirely inde- 
pendent route the same set of Maxwell's equations modified to fractals as those 
in Section 2. In analogy to (Li & Ostoja-Starzewski, 2009) which treated the 
continuum mechanics of fractal elastic media independently by Newtonian and 
Lagrangian approaches, this shows that the formulation of continuum physics 
equations based on the product measures is consistent. In the case of isotropy, 
our equations do not reduce to those of Tarasov (2010). 

5 Second order differential equations of electro- 
magnetism 

In Gaussian units Eqns (3.9), (3.3) and (3.8) are 



V • E = 0, 



cdt 



V D -B = 0, 



B + V x E = 0, 



(1) 



—sr- V xB + — J = 0, 
c ot c 

Taking derivatives gives the second-order Maxwell's equations for a fractal 



7D wDm „ D ,„ D 1 9 2 E . 4tt d 



V".(V"E)- V^(V^-E) 



or 



» 



Ek, P /c{ 



(v) 



E k ,i /c\ 



(0 



c 2 dH c 2 dt 



1 d 2 E k 4tt dJ k 

' p ~ ^~d^~ + T~dT 



(2) 



Now, consider two special cases: 

Conductor. With Ohm's law for anisotropic media (2.22), the electrody- 
namics equations simplify to 



V D • D = 0, 
V D B = 0, 

-^-H + V D xE = 0, 

c ot 



(3) 



a_ SE 



V x H = 0, 
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Upon using the fractal curl operation, we obtain a parabolic type equation 



V-.V-H=^.§ 

cr at 



or 



1 

Jp) 



H k , P /c\ 



(p) 



dH„ 
c 2 dt 



(4) 



Dielectric. From (5.3) we find a hyperbolic type equation 

cr d 2 E 



or 



Jp) 



jp l„K, 
111 kip I c i 



(P) 



d 2 t 



a km d 2 E„ 
c 2 dt 2 



(5) 



6 Electromagnetic Energy and Stress 

Pontying vector. Now, we return to (3.3) and (3.7). Multiplying the first of 
these equations by H, and the second one by E, leads to 



V D x E -H= — H 

dt 

V D x H E = J E+— D E 

dt 

Subtracting (6.2) from (6.1), on account of 

H • V D x E + E • V D x H = V D ■ (E x H), 

we obtain 

V D ■ (E x H) + J • E = --^D • E-^-B • H 

y J dt dt 

or 



1 d d 

e ijk-^j{EkHi),j + J t Ei — —Ei—Di — H l —B l 



(1) 
(2) 

(3) 



(4) 



Integrating (6.4) over the fractal's volume and applying the Green-Gauss theo- 
rem, yields 



f G • n dSd + / J • E dV D = - [ 

Js Jw Jw 



|d.e-|b.h 

dt dt 



dV D (5) 
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where 

G = E x H (6) 

is identified as the Poynting vector, which is seen to have the same form as in 
non-fractal media. 

Now, the electric and magnetic force densities are 

f E = pE + P • V D E 

(7) 

f M = J M x B + M • (B D V) = - J E x B + M • (BV D ) 

This has to be accompanied by an expression of the Ampere law in electrostatic 
units 

V ^ xH =^J + ir) = ^(J + iD) (8) 

c c c c 

7 Conclusion 

We determine the equations governing electromagnetic fields in generally anisotropic 
fractal media using the dimensional regularization together with the recently in- 
troduced product measure allowing an independent characterization of anisotropy 
in three Cartesian directions of the Euclidean space in which an anisotropic frac- 
tal is embedded. This measure also gives rise to fractal gradient, divergence and 
curl operators, which are shown to satisfy the four fundamental identities of the 
vector calculus. Next, the conservation of the electric charge on a fractal shows 
that the fractal divergence has to be used. With this Ansatz, in the first place, 
we obtain Maxwell equations modified to generally anisotropic fractals using 
two independent approaches: a conceptual one (involving generalized Faraday 
and Ampere laws), and the one directly based on a variational principle for elec- 
tromagnetic fields. In both cases the resulting equations are the same, thereby 
providing a self-consistent verification of our derivations. 

Just as in the previous works of Tarasov, we find that the presence of 
anisotropy in the fractal structure leads to a source/disturbance as a result of 
generally unequal fractal dimensions in various directions. However, our modi- 
fied Mawell equations do not coincide with those of Tarasov (2006, 2010). The 
modifications due to fractal geometry carry over to the parabolic (for conduc- 
tor) and hyperbolic (for dielectric) equations, while the Poynting vector is found 
to have the same form as in the non-fractal case. Finally, the electromagnetic 
stress tensor is reformulated for fractal systems. Overall, all the derived rela- 
tions depend explicitly on three fractal dimensions on in the respective Cartesian 
directions Xi, i = 1, 2, 3, as well as the spatial resolution; upon setting all on = 1, 
these relations reduce to conventional forms of Maxwell equations. 
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